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^ ■ Abstract 

-xf : 

This paper deals with the estimation of the performance and stability for a 
shrouded-fan unmanned rotorcraft whose mission profile also prescribes the flight 
in ground effect. The not so simple estimation of the aerodynamic coefficients and 
of the thrust in the various situations makes the performance calculation and the 

1/3 _ 

stability analysis difficult tasks. This is due to the strong interaction between the 

"1/3 

fan flow and shroud that causes quite different flow structures about the airframe 

P-i- 

depending on flight conditions. A further difficulty is related to the ground effect 
>- ■ which produces substantial modifications in the rotor thrust and aerodynamic coeffi- 

m 

cients. To evaluate performance and stability, two models have been developed. One 



determines the aerodynamic coefficients of the shroud, whereas the other one calcu- 



Q\ ', lates thrust and moment of the rotors system. Both models take into account the 

mutual interference between fan flow and fuselage and ground effect. Performance 
$_i ■ and stability are then discussed with reference to significant flight conditions. 

a ; 

Nomenclature 
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Ci, C m , C n hull aerodynamic moment coefficients 

Ci T ,Cm T ,Cn T rotor moment coefficients 

C x , C Y , C z hull aerodynamic force coefficients 

Cx t ,Cy t ,Cz t thrust force coefficients 

Ct = —Cz Tl — Cz Tl rotors system thrust coefficient 

Cq = C nTl — C nTu rotors system torque coefficient 

d rotor diameter 

D overall shroud diameter 

g acceleration due to gravity 

h rotorcraf altitude 

I x , I y , I z principal moments of inertia 

I xz product of inertia 

J merit function 

V i • 

k = velocity parameter 

w D 

L, M, N hull aerodynamic moment components 

Lt, Mt, Nt rotor moment components 

Q = Nti — Nt u rotors system torque 

R rotor radius 
tcD 2 

S = rotorcraft reference surface 

4 

u = (5a, 8b, Sc, dp) control vector 

v = (u,v,w) inertial velocity in body axes 

v local fluid velocity 
v 

v = — = (u, v, w) dimensionless velocity in body axes 
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V 

wd = Wi - V sin a 
W 

Wi 

x,y,z 

xb,Vb, z b 
X,Y,Z 
X T , Y T , Z T 
x 

Greek symbols 

a,(3 x 

a x ,P 

7T 

a = - - a x 

5a 
5b 
5c 
5 P 

e 

^ = (A x , \y, A 2 ) 



velocity modulus 

fluid velocity through the disk 

weight 

induced velocity by the rotor flow on the rotor disk 

inertial coordinate 

coordinate in body axes 

aerodynamic force components in body axes 

thrust force components in body axes 

position 



aerodynamic angles 

angles of attack and sideslip, respectively 

axial angle of attack 

lateral cyclic control, positive right 

longitudinal cyclic control, positive after 

collective control, positive up 

differential collective control, positive right 

ground inclination with respect to the horizontal plane 

Euler angles 

ground normal unit vector in body axes 



advance ratio 



v w 






dimensionless induced velocity at the disk 
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required power 
p air density 

r = yjl + A(h/R) 2 — 2h/R ground effect perturbative parameter 
(p, q, r) angular velocity vector in body axes 

Q rotors angular velocity 

Subscript 
a.c. aerodynamic center 

e.g. vehicle center of gravity 

u upper rotor 

1 lower rotor 

oo value calculated out of ground effect 

Introduction 

Performance and stability calculation of a shrouded-fan unmanned aerial vehi- 
cle (UAV) is a very difficult task for various reasons. One of these is related to 
the intensive aerodynamic interaction between shroud and rotors 1,2 that, because 
of the large excursion of angle of attack and fans system induced velocity, deter- 
mines quite different flow fields about the UAV 3-9 . The simultaneous variations of 
angle of attack, flight speed and the rotors induced velocity, determine a wide set 
of flight conditions, where the aerodynamic and propulsive force and moment coef- 
ficients exhibit large variations. Furthermore the vicinity of the ground makes this 
study more complex, due to the influence that the ground exerts on both fuselage 
and rotors system. Therefore, the performance and stability analysis requires an 
accurate vehicle aerodynamics that calculates all the aerodynamic coefficients in a 
wide range of angle of attack and fans working regime, and an adequate model to 
determine the rotors force and moment in the different flight conditions. The UAV 
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here considered is a shrouded-fan made of a toroidal airframe at the center of which 
are placed two counter-rotating rotors driven by three two-stroke air cooled engines 
(see Fig. 1), the main characteristics of which are reported in Tab. I. This vehicle 
is similar in shape to that studied in Ref. 1 and is the results of a project jointly 
accomplished from the University of Rome "La Sapienza" 2 and the Polytechnic of 
Turin 1 . The fuselage geometry is the result of a parametric study which analyzed the 
lift capability of different toroidal geometries in the presence of the rotor flow. The 
selected fuselage cross section develops the highest lift in hovering for an assigned 
hull volume. The main function of the hull consists in to incorporate the avionics, 
fuel, payload and the possible flight related hardware. From an aerodynamic per- 
spective, the shroud should contain the streamlines to follow the duct and exit at an 
adequate velocity imposed by the exit area. Moreover the fuselage causes an intense 
suction force, which is produced by the inlet lip, that makes the fuselage a body with 
a nonnegligible lift capability. Hence the shroud generates a sizable improvement of 
the static thrust which depends on the reciprocal interaction between fans system 
and fuselage 5 . During axial flights, the advantage in static thrust furnished by the 
shroud rapidly diminishes 5 as soon as the UAV increases its speed. 

For what concerns the ground proximity, it is responsible for several effects. One 
of these is represented by the fuselage aerodynamic coefficients variations which 
occur in hovering. In such situation ip — $ — 0, then the ground imposes a bound- 
ary condition on the aerodynamic field that symmetrically reduces the downwash 
of the rotors. The consequence is that the aerodynamic coefficients depend on the 
dimensionless height h/R, and an important loss of the fuselage lift capability is 
observed 10,11 . If the vehicle flies at an arbitrary inclination with respect to the 
ground, the interaction between rotors flow and ground alters the rotorcraft aero- 
dynamics in such a way that the pressure distribution on the airframe depends on 
the height and on the vehicle attitude. Hence the aerodynamic coefficients, which 
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are functions of the aerodynamic angles, in ground proximity also depend on the 
dimensionless height and on the Euler angles. 

The vicinity of the ground also influences the rotors system. It has been observed 
that for an unducted rotor in hovering the thrust considerably increases as the height 
diminishes 4-9 . Furthermore the attitude modifies rotor actions 12 so that, in ground 
effect thrust and torque depend on the Euler angles. As a result, the attitude can 
play an important role for the determination of the flying qualities depending on 
the height. Also the ground inclination with respect to the horizon is of great im- 
portance in ground effect. This can be the case of a rotorcraft that approaches a 
ship flight deck 12,13 whose inclination influences the vehicle aerodynamics and ro- 
tors thrust and moment 12 . The deck oscillations due to the ship attitude motion 
determine unsteady thrust and aerodynamic force and moment which influence the 
rotorcraft dynamics. Another source of the significant variations of the performance 
and stability in ground effect is nonzero flight path angle which determines incre- 
mental aerodynamic force and thrust 14 . This may happen during maneuvers that 
are used to avoid collisions with low obstacles. The aforementioned effects are more 
pronounced at low height where the ground effect is more intense and can determine 
sizable variations of trim controls and rotorcraft stability margin. 

Although several studies on the mutual interference between hull aerodynam- 
ics, rotor flow and ground effect 12 ' 14-17 were developed, to the author's knowledge 
a general analysis of these phenomena in terms of rotorcraft flight dynamics has 
not received due attention. Therefore, the objective of the present work is to de- 
velop an accurate mathematical model that may be adopted for the analysis of the 
performance and flying qualities of the UAV. 

In the present study two models are proposed. The first one calculates the shroud 
aerodynamic coefficients, whereas the second one determines the rotors force and 
moment through the blade-element theory. Both models account for the mutual 
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interaction between hull and rotors and the ground proximity. The aerodynamic 
model exhibits free parameters that have to be identified. These models allow the 
forces and moments to be accurately formulated and the performance and flying 
quality to be interpreted. Obtained results are compared with data reported in 
the literature, and finally a detailed analysis of the performance characteristics and 
flying qualities of the vehicle is carried out. 

Aerodynamic Model 

A method to calculate the aerodynamic force and moment on the UAV shroud 
configuration, is now presented. This procedure, based on physical considerations, 
takes into account the reciprocal influence between rotors flow and fuselage and the 
possible presence of the ground. First, it is necessary to establish the flow structure 
about the airframe that in turn directly depends on the rotors system working 
regime. Fig. 2 shows the regions of the four characteristic working regimes in terms 
of the induced and axial velocities 3-6 . In this study the vehicle aerodynamics is 
modelled according to the normal working state which corresponds to situations 
where the induced velocity is sufficientely high with respect to the axial velocity 
V^sina. In this situation the rotors flow enters from the top of the rotorcraft, the 
air is then ejected towards the bottom and the power is transfered from the rotors to 
the air. The other working regimes, such as the windmill-brake state, the turbulent 
wake state and the vortex ring state are not studied in the present analysis since 
they are considered to be off-design situations which correspond to more complicate 
flow conditions that strongly reduce the hull lift capability. 

Figure 3 shows the rotorcraft in normal working state, out of ground effect. The 
rotors flow and the external stream are separated by the wake W which originates 
from the separation line T that is assumed to be equal to the circumference of the 
exit area. Because of the axial symmetry around zb, the sideslip angle does not 
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influence the shroud aerodynamics. Aerodynamic force and moment lie on 
the plane zb — V which is defined by z B and v. This plane defines a reference frame 
(x v , y v , z v ) whose axes are oriented so as z v = zb, y v is normal to both z v and 
v, whereas x v is perpendicular to both y v and z v . In this frame and are 
expressed as 

F A = (X,0, Z) 

(1) 

Qa = (0, M, 0) 

where all the components are independent on the sideslip angle. Due to the axial 
symmetry, instead of using the angle of attack a x and sideslip f3, it is convenient to 
introduce different aerodynamic angles, i.e. a and f3 x , that are defined as follows 1 ' 2 



u 




u 




cos a cos (3 X 


V 


= V 


V 


= V 


cos a sin f3 x 


w 




w 




sin a 



The components in body axes of F^ and in terms of (3 X are expressed as 
X = X cos j3 x , Y = X sin (3 X , Z — Z, 

(3) 

L = -M sin /3 X , M = M cos f3 x , N = 
The hull aerodynamics is strongly influenced by the rotors flow. In this model 
an equivalent actuator disk is placed in the mid of the two rotors, whose induced 
velocity Wi coincides with that produced by the rotors at the actuator disk plane. 
Therefore, the rotor flow is supposed to be generated by an uniform layer of doublets 
with a proper value of Wi. The induced velocity, that in turn represents the effect 
of the collective pitch on the vehicle aerodynamics, is determined by means of the 
rotor model described later in the dedicated Section. Non-uniform distributions of 
Wi on the disk plane, caused by the cyclic pitch, are neglected. 

To derive the analytical expressions of the aerodynamic force and moment coef- 
ficients, the rotorcraft is considered in a steady state potential flow. Then, the flow 



8 



velocity in a point x of the space is given by 18 

Ov 3d 

V ^ ) = d^ V+ d^ WD (4) 
where dv/dv and dv/dwo are influence matrices. In Eq.(4) the first term gives 
the effect of the flight velocity on the vehicle aerodynamics, whereas the second one 
provides the contribution caused by the rotors flow. Force and moment depend on 
the local pressure distribution which is related to the square of the local velocity 

v 2 = (hi cos 2 a + H2 sin 2a + /X3 sin 2 a)V 2 + 2Vwd(xi cos a + X2 sin a) + g 2 wo 2 (5) 

by means of the Bernoulli theorem, where fa and Xi are functions of the position 
through dv/dv and dv/dwo- Because of the symmetry around zb, the pressure 
does not depend on the sideslip angle. The aerodynamic force and moment are then 
expressed as surface integrals of the local pressure 

Fa = I Is Pn dS ' Qa = I U T ~ XnpdS (6) 

Substituting Eq. (5) into Eqs. (6) one obtains the aerodynamic coefficients on the 

plane z v — V in terms of k and a 

Cx{oc, k) = Fi(k, a)Xi sin 2a + F 2 (k, a)X 2 cosa 

C z (a, k) = F (k, a)Z + F 1 (k)Z l sin a + F 2 (k, a)(Z 21 cos 2 a + Z 22 sin 2 a) (7) 

C m (a, k) = Fi(k, a) Mi sin 2a + F 2 (k, a)M 2 cos a; 
that are defined by means of 



F A = ^pV 2 S (CxA C z ) 
Qa = \pVfSD (0, C m , 0) 



(8) 



where 



Vf = \ju 2 + v 2 + (w — Wi) 2 = \Jv 2 cos 2 a + w 2 D , (9) 
while the functions Fi, F 2 and F 3 are given by 

F o(M) = T2 \~— T , Fi(k,a) = — , F 2 (k,a) = \ (10) 

k cos^ a + 1 k l cos^ a + 1 k 2 cos z a + 1 
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The aerodynamic coefficients in body axes are then obtained as the projection of 
Eqs. (7) in body frame, i.e. 

C x = C x cos(3 x , C Y = C x sm(3 x , C z = C z , 

(11) 

C\ = —C m sin (3 X , C m = C m cos(3 x , C n = 
The angular velocity is taken into account by rotational derivatives, the expressions 
of which are 18 

Cx p = 0, Cx q = M 33 sin a, Cx r — —Mu cos a sin (5 X 

C Yp = — M 33 sin a, C Yq = 0, C Yr = M n cos a cos (3 X 

C Zp = Mu cos a sin f3 X) C Zq = — M n cos a cos (3 X) C Zr = 

Ci p = K12 sin a, 

Q q = (Kas-Ku) sin a, (12) 
Ci r = (K 33 — Ku) cos a sin f3 x — K\ 2 cos a cos f3 x 
C mp = (Ku- K 33 ) sin a 
C mq = K12 sin a, 

C mr = (Ku — K 33 ) cos a cos j3 x — K w cos a sin f3 x 

where 

Mij = Ma F (k, a) + Mij ^(k, a) 

(13) 

K i:j = K i:j Q F (k, a) + Kij iFi(k, a) 
Eqs. (12) represent the expressions of the rotational derivatives of a body which 
exhibits an axial symmetry around zb- The quantities Mj, X iy Zi appearing in 
Eqs. (7) and Z^, , 1, and 1 appearing in Eqs. (12- 13) are free 
parameters which depend on the vehicle geometry. 

In ground effect, the rotor downwash is limited by the ground surface, so that 
the velocity field about the airframe is modified with respect to Eq. (4). To study 
this situation consider now Fig. 4, where the rotorcraft in close proximity to the flat 
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ground is analyzed. The wake W is considered to be a rigid cylinder whose length is 
assigned in advance in such a way that it does not intersect the ground plane. The 
ground inclination with respect to the rotorcraft is represented by the unit vector 

A = (A-r, \ y , A 2 ) = (— sini?, simp cos-#, cos (p cos i?) (14) 

which is the normal unit vector of the ground plane in body axes. The presence of 
the ground is simulated by placing a specular image of the rotorcraft at an equal 
distance below the ground plane. According to the potential flow theory 18 , the local 
velocity is the sum of those induced by the vehicle and by its mirror image, i.e. 

dv d 2 v 

" (x) = ^ (X) + W^} WdX + dXd^ Xv (15) 

The first addend is the local velocity induced by the entire rotorcraft, whereas the 
second and the third ones provide, respectively, the interaction between rotors flow 
and ground plane and the influence that the ground exerts on the fuselage. The 
pressure about the airframe is then altered by the vicinity of the ground, so that 
the aerodynamic coefficients also depend on height and attitude. Now, to obtain 
the expressions of the force and moment in ground effect, the kinetic energy of the 
stream T is considered. Following Ref. 18, T is a definite positive quadratic form 
that is a function of v, wc and A. In order to determine the expression of T, consider 
now two sets of flight conditions. The first of these corresponds to horizontal flights 
[ A = (0, 0, 1) ], where v describes a circular cone about z B , whereas the second one 
consists of axial flights [ v = (0, 0, V) ] with A which makes a circular cone about 
z B . In both cases, due to the axial symmetry about z B , the kinetic energy must 
remain unaltered. As a result the kinetic energy of the stream can be written in the 
form 

T = f^u 2 + v 2 ) + f 2 w 2 + h w D (u\ x + v\y) + U w D w\ z (16) 
where /i, fa, fs and f± are functions of h/R that, according to Ref. 19, are written 
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in terms of h/R by means of the perturbative parameter r, i.e. 

fk=j2f km r m , (k = 1,2,3,4) (17) 

m=l 

The shroud aerodynamic force and moment increments developed in ground effect 
are then expressed using the Lagrange equations method 18 

d dT d dT 
AF = (AX, AY,AZ) = ------V^ 

(18) 

dT 

AQ = (AL, AM, AN) = -v x — + r AG x AF 

ov 

where r AG = r a.c. _ r c.g. = ( x ag 5 Yag; 0) is the position of the aerodynamic center 
with respect to e.g., which varies in ground effect. Out of ground effect it is assumed 
that r a c . = r c . 9 . while close to the ground, r a . c . depends upon the flight conditions. 
Since the yaw moment must be equal to zero in all the situations, the projection 
along zb of r AG x AF must be identically equal to zero. The consequence is that 
r AG assumes the form 

^ = C(^,^,0) (19) 

f< ^ 'oo "^OO 

where ( is a quantity depending on h/R that is expressed through r 

C = E X k r k (20) 

k=0 

Substituting Eq. (16) into Eqs.(18) one obtains the incremental force and moment 
in ground effect 

AX = -(2/xm + f 3 w D X x ) - (2f lh u + f 3h w D \ x ) Vsin-f 

AY = -(2f lV + hw D \ y ) - (2f lh v + h h w D \ y ) V sin 7 (21) 

AZ = -(2f 2 w + Uw D \ z ) - (2f 2h w + U h w D \ z ) Vsin'j 

AL = -2(/ 2 - fi) vw - f 3 w D X z v - xR[{2fiv + f 3 w D X y ) + (2f lfl v + f 3h w D X y ) V sin-/} 
AM = 2(/ 2 - h) uw + f 3 w D X z u + X R[(2fiU + /a^A,) + (2f lh u + /s^dA.) ^sin 7 ] (22) 
AN = 
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where the quantities f\ and fh are 

fk = fkvV + fk WD w D , f kh = fwV + fkw D w D (23) 

Following the classical theory of the lifting bodies, both f k and f h are proportional 
to the circulation around the hull. In particular, f k is the time derivative of 
which is caused by the wake, while fkh is the derivative of f k with respect to h 
which is related to the sink rate. Accounting for the variations of ( and fk, the 
aerodynamic force coefficients can be written in the form 

2 3 

AC X = Yl F m{k, a)r n [(A mn u + B rnn \ x ) + (C mn u + D mn \ x ) sin 7] 

m=0 n=l 
2 3 

AC V = E E F m(K a)r n [(A mn v + B mn \ y ) + (C mn v + D mn \ y ) sin 7] (24) 

m=0 n=l 
2 3 

AC Z = E E F ™(k, a)r n [{E rnn w + G mn A z ) + (H mn w + L mn \ z ) sin 7] 

m=0 n=l 

whereas the aerodynamic moment coefficients are 

2 3 

ACi = E E ^m(fc, «)r"[(iV m „{) + P mn {)w + Q m „X z v + R mn X y )+ 

m=0 n=l 

(SmnV + T mn wu; + U mn X z v + V mn \ y ) sin 7] 

2 3 ^5) 
AC m = - E E F m (k, a)r n [(N rnn u + P mn uw + Q m „X z u + R m „X x )+ 

m=0 n=l 

) sin 7] 

AC n = 

In Eqs. (24) and (25) the coefficients A mn , B mn , C mn , D mn , E mn , G mn , H mn , 
Lmni N mn , P m n) Qmni Rmni ^mni F mni U mn and V mn are the free parameters of the 
aerodynamic model in ground effect. Since for k — > 0, AC X = AC y = AC t = AC m 
do not depend on the aerodynamic angles, one has A 0n = Co n = E 0n = H 0n = N Qn 
= Pon = Qon = Son = T 0n = U 0n =0. The terms with sin 7 yield the influence of the 
sink rate on the shroud aerodynamics. A non zero sink rate determines an unsteady 
ground effect that modifies the aerodynamic coefficients according to the Bernoulli 
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theorem 

+ ~^r- + - = const 26 
at 2 p 

where 86 /8t = 86 /8h h = d(f)/dh V r siri7 is the unsteady term that produces the 
dynamic ground effect. As for the influence of the spin rate, Eqs.(12) are considered 
to be valid also in ground effect. 

Rotors System Model 

The rotors actions provide lift force and control moments to manage rotorcraft 
attitude. Pitch and roll are controlled through longitudinal 5 b and lateral 5 a varia- 
tions of blade pitch, whereas the yaw control is carried out by means of differential 
variation 5p of the collective pitch on both rotors, whose angular velocity is kept con- 
stant by a RPM governor. The blade pitch is controlled by a mechanism consisting 
of two independent swash-plates, each driven by three actuators. 

Thrust forces and moments developed by the rotors are given by the equations 6 
(X T , Y T , Z t ) =-n P Q 2 R 4 (C XT , C yT , C ZT ) 

(27) 

(L T ,M T ,N T ) = npQ 2 R 5 (Q T ,C mT ,C nT ) 
where C XT ,C yT ,C ZT and Ci T ,C mT ,C riT are, respectively, thrust and torque coeffi- 
cients. They are obtained with the blade-element theory which calculates thrust 
and moment through analytical integration of the aerodynamic load along the blade 
span assuming steady-state aerodynamics 1 , where the effects of the blade-tip losses 
and the mutual influence between the two rotors are neglected. The determination 
of the rotors regime is made by imposing that the thrust coefficient C ZT (7 X , 7^, 7 Z ), 
calculated with the blade-element theory, is equal to that obtained by the actuator 
disk theory 20 , that is 

C ZTa = 2Ai/y2 + A?-2 72 A, sin a (28) 
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where 



and 



_ U + U s _ V + V S _ W + W s Wi_ 



. . dv 

{u s ,v s ,w s ) = — ■ v (30) 

is the velocity induced by the shroud at the center of the actuator disk plane. As a 
result, one obtains the following equation 

C«r(7*,7»»7*) = + \ 2 - 2 lz ~K sin a ( 31 ) 

the solution of which gives the dimensionless induced velocity A, and the corre- 
sponding thrust coefficient. To take into account the ground effect, the ground is 
modelled by means of the mirror image of the entire rotorcraft. Equation (31) is 
still considered valid even though (u s ,v s ,w s ) now accounts for the induced velocity 
of the mirror image calculated in the origin of the actuator disk plane, i.e. 

, . dv dv d 2 v 

(M - V - Ws) = d^- V+ 8X^ WdX + dXd^ Xv (32) 

A nonzero sink rate determines an unsteady ground effect for the rotors system, 
then following the Bernoulli theorem, the pressure distribution along each blade is 
augmented by the quantity 

Ap= -p^Vsin-f (33) 
that expresses the pressure increment caused by the dynamic ground effect. 

The equations of motion 
A full nonlinear six-degree-of-freedom model of the vehicle is now defined. Fol- 
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lowing Etkin 21 , the equations of motion for the rotorcraft are written as follows 

u = q — — q sin v — q w + r v 

W 

■ Yt + Y i 9 • 

v = g — — h g cos v sm ip — r u + p w 

Zt + Z 

w = g — — h g cos v cos ip — p v + q u 

(p = p + sin ip tan •& q + cos ip tan i? r 

i? = cos ip q — sin ip r 

ip = sin ip sec i? g + cos <p sec $ r 



[4*P<7 + (4 - h)qr + L + L T ]I Z + [-I xz qr + (4 - + iV + iVr]/^ ^ 



I I - I 2 

1 X 1 Z J K 



. hz {r 2 - P 2 ) + (I z -I x )pr + M + M T 
Q = J 

. = [-I xz qr + {h ~ I y )pq + N + N T }I X + [I xz pq + (I y - I z )qr + L + L T }I XZ 
T I I - I 2 

1 x 1 z 1 XZ 

x = (cos d cos ip) u + (sin 93 sin $ cos ^ — cos ip sin ■0) t> 
+ (cos ip sin $ cos ip + sin sin ip) w 

y = (cos sin ijj) u + (sin ip sin sin -0 + cos p cos -0) v 
+ (cos (/? sin sin ^ — sin ip cos -0) w 

i = (— sin-i?) -u + (sin(pcosi9) v + (cos </? cos u> 

It is worth to remark that both the coaxial rotors cause moments that are transmit- 
ted to the fuselage. Since the rigid rotors have the same moments of inertia with 
respect to the rotation axis, their gyroscopic effects are balanced each other and 
therefore do not appear in the rigid body moment equations. 

Identification of the Aerodynamic Model Parameters 

This section describes the identification procedure for the free parameters of the 
aerodynamic model. The algorithm consists of an optimization method, based on 
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the least-square procedure, that is represented by the problem 

J = Y.[( C k)cFD-C k ] 2 = mm (35) 

k 

where is the generic aerodynamic coefficient calculated with the proposed model, 
whereas {C^)cfd is the same coefficient computed by CFD simulations. These sim- 
ulations are obtained using VSAERO by Analytical Methods, Inc. 22 , which is a code 
based on a boundary integral formulation. The code accounts for the aforementioned 
effects such as the interference between fan flow and airframe, the ground effect and 
the presence of the separated flows enclosed by the wake. 

The design variables of the problem given by Eq. (35) are defined as the argu- 
ments of J, that are the aerodynamic model free parameters. 

The free parameters of the model out of ground effect are determined first. They 
are the coefficients of the Eqs. (7) and (13) which are identified through several 
CFD calculations. Once the simulations are carried out, the quantities dv/dv and 
dv/dwD are also determined. The results are shown in Fig. 5, that reports drag, lift 
and pitch moment coefficients, where the continuous lines and the solid symbols are, 
respectively, the results calculated by the model and the data obtained by VSAERO. 

For what concerns the ground effect, the free parameters are the coefficients of 
Eqs. (24) and (25). To assess the simultaneous influence of a, f3 X) (p, h/R and k 
on the vehicle aerodynamics in ground effect, a large number of CFD simulations is 
made with VSAERO, for several flight conditions. The model free parameters are 
then calculated using the aforementioned minimization algorithm, so that also the 
influence functions dv/dXwo-, d 2 v/d\dw and d<p/dh are evaluated. 

Validation of the Models 

To validate the proposed models, the present results are compared with some 
existing data in the literature. 
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In Fig. 6 the aerodynamic coefficients in body axes in terms of a are depicted. 
The present results (continuous lines) are compared with those from Ref. 23 (solid 
symbols). This latter model is based upon Fourier expansion of the vorticity distri- 
bution that accounts for the effects of camber, taper and thickness. These data are 
determined by combining all these effects to obtain a similar toroidal fuselage to the 
shroud here studied. Some discrepancy is apparent, especially for C z when k tends 
to zero, that could be caused by the difference between the two geometries. In spite 
of that, the present results can be considered in somewhat comparable with those 
obtained using the model of Ref. 23. 

As for the ground effect, it is known that for a platform in the presence of a lifting 

jet 10 , the hull exhibits a substantial lift reduction 10,11 that becomes more influent 

as the fuselage approaches to the ground. In the case of hovering (V = ip — d 

=0), it can be shown through the Buckingham theorem that, for given thrust and 

hull geometry, the dimensionless parameters that describe the lift loss phenomenon 
AZ 

are — — and h/(D — d) 10,n . Fig. 7a shows the present results (continuous line) 

Zoo 

in comparison with the data from Ref. 10 and 11 (dashed line and solid symbols, 
respectively), that are static measurements of the lift losses in ground effect for 
a circular platform. It is worth to remark that the present data are referred to 
a toroidal fuselage, while Refs. 10 and 11 deals with circular planforms with a 
centrally-located lifting jet. Nevertheless, the obtained results are qualitatively in 
good agreement with the those of Refs. 10 and 11. 

The aerodynamic coefficients in ground effect also depend on the vehicle attitude 
and sink rate. Then, Fig. 7b, 7c and 7d show the increments of the longitudinal 
aerodynamic coefficients (continuous lines), calculated at ip — 0, k — for various 
h/R, while Fig. 7e gives the derivative dC z /d^j at a = rc/2, v.s. k, where each curve 
represents a given h/R. The results are represented together with those of a rigid 
vortex ring in the vicinity of a wall (solid symbols), where the vortex geometry is 
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assumed to be equal to the circumference defined by the centers of the shroud cross 
sections, whereas the vortex circulation is selected in such a way that the induced 
velocity at the center of the vortex ring is equal to Wi. As well as in the cases of 
the hull and rotors system, the wall is simulated by means of the vortex ring mirror 
image, where the induced velocity on each vortex element is calculated by means 
of the Biot Savard law, whereas the corresponding action is determined using the 
Kutta-Joukowsky theorem. Some disagreement is evident, due to the geometrical 
differences between toroidal fuselage and vortex ring, that becomes more significant 
as h/R goes to zero. Nevertheless, the obtained results seem to be in adequate 
agreement with those of the vortex ring. Because of the symmetry around zb, AC x 
and AC m are odd functions, while AC Z is an even function of AC Z varies in 
accordance to Fig. 7c and presents minor variations with t? but for h/R — > 0, while 
C m exhibits a negative slope whose absolute value increases as h/R approaches to 
zero. This last characteristic corresponds to a positive contribution to the rotorcraft 
stability in ground effect. The plot of Fig. 7e, which illustrates the effect of the 
sink rate on the hull aerodynamics, shows the derivative dC z /d^ in terms of k. A 
negative flight path angle (positive sink rate) produces a negative variation of C z 
that corresponds to a an increment of the shroud lift coefficient. For each curve, 
dCz/d'-) exhibits its maximum value at about k=l. As for the lateral coefficients, 
because of the hull symmetry, when t? is changed with ip, then C x , C m are changed 
in C y and -Ci, respectively. 

For what concerns the rotors model, various calculations have been carried out. 
Fig. 8 gives Ct and Cq of the rotors system in axial flight with (continuous lines) 
and without (dashed lines) shroud in terms of the advance ratio for several collective 
pitch. In both cases, the rotors system exhibits the behavior of a propeller in axial 
flight, whereas the presence of the shroud produces sizable variations in the thrust. 
For each collective pitch, at 7 = (static case), the fuselage generates a thrust gain 



19 



of about 25% with respect to the unshrouded rotor, whereas, as soon as the advance 
ratio increases, a more limited thrust gain is observed. As for the torque coefficient, 
minor variations are observed. This results are in accordance with Refs. 3 and 5, 
where similar configurations of ducted propeller are treated. 

The vicinity of the ground significantly modifies the rotors characteristics. The 
plots in Fig. 9 show the thrust coefficient, induced velocity and induced torque coef- 
ficient v.s. h/R. In Fig. 9a the solid symbols are from Ref. 4, which represent flight 
tests accomplished with different helicopters, whereas the continuous line yields the 
present results that are obtained by applying the proposed rotor model to a single 
free rotor. Fig. 9a also shows the Cheesman and Bennett results 16 (dot-dashed line), 
wherein the presence of the ground is modelled by placing under the ground a mirror 
image which consists of a simple fluid source whose mass flow coincides with that of 
the rotor. The dashed lines give the data obtained using the Hayden 17 method which 
estimates the influence of the ground in hovering through flight test measurements. 
According to the literature 4 ' 5 , the Hayden results are found to overpredict the rotor 
thrust. The figure shows that the present results are in good agreement with the 
aforementioned flight tests measurements. Fig. 9b and 9c compare the induced 
velocity and the torque coefficient with the other data sources. The solid symbols in 
Fig. 9b are from Ref. 12, that represent the normalized average inflow in terms of 
h/R. According to Xin and Prasad 12 , the continuous line is calculated by applying 
the proposed model to rotor of the Yamaha R-50, which is a small size, remotely 
piloted helicopter with a rotor diameter and a rotor speed of 3.07 m and 850 R.P.M., 
respectively. The present results, which are compared with Ref. 12 (solid symbols, 
see figure), shows that the rotor model gives results in good agreement with those 
reported in the literature. The continuous line with the filled symbols provides the 
data for the UAV here studied. 

Also the attitude influences the rotor forces and moments in ground effect. In the 
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diagrams of Figs. 10a, 10b and 10c, the coefficients C Zt , C Xt and C mT at V=0 are 
reported as functions of where continuous an dashed lines represent the data for 
shrouded and free rotors, respectively. As well as in the case of the hull, because of 
the vehicle axial symmetry, Cz T is an even function, while Cx T and C niT are both 
odd functions of the pitch angle. While the pitch angle produces minor variations 
on Cz T , C mT exhibits negative slopes whose absolute value increases as h/R tends 
to zero. This last characteristic is the contribution of the fans system to the vehicle 
stability in ground effect. The presence of the shroud significantly modifies the 
thrust force and moment coefficients. In particular a sizable Cz T increment occurs, 
which in turn is in agreement with Fig. 8, and an important slope reduction of both 
Cx T and C mT is observed that is caused by the shroud induced velocity. 

Again, following Ref. 12, the present model is now applied to the Yamaha R-50 
rotor (Fig. 10 d). The curves give the torque coefficient vs. h/R at different 
where the solid symbols are from Ref. 12, while the continuous lines correspond to 
the present results. According to the literature 12 , for each h/R, the vehicle attitude 
tends to reduce the torque coefficient, and the comparison with Ref. 12 shows that 
the maximum difference between the two methods results to be always less than 7%. 

Finally, to assess the dynamic ground effect of the rotors system, the thrust vari- 
ations in term of the dimensionless sink rate are shown in Fig. lOe. Continuous 
and dashed lines give, respectively, the present model and the results obtained with 
the unsteady actuator disk theory 20 , whereas the symbols represent the thrust cal- 
culated by a code based on the vortex lattice method 24 , that simulates the dynamic 
ground effect by means of the mirror image of the actuator disk that moves with a 
flight path angle equal to -7. These comparison demonstrates the good agreement 
between the present data and the results obtained by the actuator disk theory and 
by the CFD simulations. 

Results and Discussion 
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In this section, to analyze the rotorcraft performance and stability in and out 
of ground effect, some significant situations which correspond to straight and level 
flight at the trim are studied. Eqs. (34) are used as equations of motion and 
the contribution of the control variables to vehicle dynamics is taken into account 
through the rotors model. The trim calculation is made by solving a minimization 
problem with assigned constraints 25 , whereas the vehicle stability is investigated by 
means of the eigenvalues analysis applied to the linearized motion equations. 

Fig. 11 depicts some of the significant variables calculated at trim. In Fig. 11a 
lift and drag coefficients are shown in terms of forward speed. Due to the hull lift 
cabability caused by the suction force developed by the inlet lip, a sizable nonzero 
lift coefficient is observed at low velocities, which is about 25 % greater than that of 
the UAV analyzed in Ref. 1. The diagram in Fig. lib shows the hull trim moment 
as the function of the flight speed. This is important for the longitudinal stability 
and controllability, since the aerodynamic moment at the trim is balanced by the 
rotors control moment, generated by longitudinal cyclic pitch. The present results 
calculated out of ground effect (dashed lines) can be compared with the data from 
Refs. 1 and 26. In Ref. 1, where a quite similar rotorcraft is studied, the trim 
moment is a rising function of the forward speed in the speed range -j- 30 m s _1 , 
and the moments calculated for the velocities of 10 and 20 m s _1 , are about 100 and 
180 N m, respectively. Ref. 26 deals with the Cypher, an uninhabited rotorcraft 
developed by Sikorsky Aircraft Corporation which is also made by a toroidal fuselage 
with at the center two coaxial rotors. There, the aerodynamic moment rises until 
to a speed of about 13 m s" 1 , where exhibits its maximum value of about 200 N m. 
Although the two vehicles treated in Refs. 1 and 26 exhibit differences with respect 
to the present one, the corresponding data seem to be in somewhat comparable with 
the results shown in Fig. lib. 

As seen, the ground proximity causes significant variations on the vehicle aero- 
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dynamics and on the fans regime in such a way that both aerodynamic and thrust 
coefficients also vary with height and attitude. This effect results to be more pro- 
nounced at lower h/R. According to the lift reduction for platforms with lifting jet 10 , 
the fuselage lift capability diminishes as the rotorcraft approaches to the ground. 
The plots in Figs. 11c, lid, lie and 1 If show the controls expressed in radians in 
terms of forward speed. For each h/R, the collective pitch 5 C (Fig. 11c) presents a 
low speed region where it remains practically constant until to the speed of about 20 
m s" 1 , over which it becomes a rising function of V, whereas for what concerns the 
longitudinal cyclic pitch 5b (Fig. lid), its values are directly related to the moment 
calculated at the trim. For h/R— 2, 3, 5c varies in accordance to the rotors thrust 
increments, whereas 5b diminishes until to 20 m s" 1 , where reaches its minimum. 
At higher velocities, according to the hull moment developed at the trim, 5b rises 
with V. The calculation made at h/R — 1 shows that 5 B is a decreasing function of 
V in the entire speed range, and for what concerns 5a and 5p, quite small variations 
are obtained, but for h/R < 1. In this last situation the simultaneous effects of 
pitch angle and height generate a strong flowfield perturbation on the two rotors 
in different fashions which must be balanced by the differential collective pitch. At 
h/R = 1 the peripheral parts of the airframe can touch the ground for certain pitch 
angles, thus the trim results to be feasible in a more limited speed range. Therefore 
the corresponding curves are broken at a flight speed less than 24 m s" 1 , where d ~ 
35°. 

For each h/R, the law i?(V) directly depends on the forces calculated at the trim 
(see Fig. 11 g). Because of the smaller variations of lift and drag with respect to 
the thrust, the pitch angle exhibits relatively small variations as h/R changes. The 
required power Il n and the velocity parameter k are represented in Fig. llh in terms 
of flight speed. Some comparisons can be made out of ground effect (dashed lines) 
with the data from Ref. 1. Close the hovering, Tl n is rather similar to that calculated 
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in Ref. 1, whereas for V ^ some discrepancies are apparent which are caused by 
the different fuselage geometries that develop diverse aerodynamic actions. The 
higher shroud lift capability than that of the vehicle dealt in Ref. 1, yields smaller 
required power when V ^ 0. The ground effect causes important modifications on 
Il n and k. At h/R = 1, for V < 10 m s" 1 , a reduction of the 10 % of IT„ is observed, 
while at higher velocities, where a d ^ 0, the interaction between the rotorcraft 
and ground generates seeable variations on both aerodynamic and thrust coefficients 
which in turn are responsible for the IT n variations for V > 10 m s -1 . The law k(V) 
calculated out of ground effect (dashed line), is similar in shape to that of Ref. 1 
where k is about 10 % greater than that here calculated. The downwash reduction 
due to the ground effect, causes smaller values of k as h/R approaches to zero. 

Next, to asses the influence of the ground inclination on the trim, the calculations 
are repeated in the cases where the ground wall is inclined with respect to the 
horizon. The wall inclination, which is here obtained by rotating the ground surface 
of an angle e around the x axis, is represented by the wall normal unit vector, that 
in the inert ial frame is given by 

A = (0, sine, cose) (36) 

To avoid possible touchdown, the trim, whose results are summarized in Fig. 12, 
is made at h/R=l for a ground inclination of e — ± 20°. In such situation the 
peripheral parts of the airframe are very close to the wall, thus the interaction 
between rotorcraft and ground determines high rotors and fuselage moments that 
yield not feasible trim for V< 20 m s -1 . Although the state variables result to be 
rather similar than those just depicted in Fig. 11, the controls are quite different 
with respect to those obtained for e =0. Relevant variations associated to 5b, 5a 
and 5p can be observed, especially for what concerns 5 a- 

The rotorcraft stability is next evaluated through the eigenvalue analysis applied 
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to the linearized motion equations. The root locus is calculated for various h/R in 
level flight varying the velocity from to 35 m s _1 (see Fig. 12 abed), whereas the 
several modes are recognized by means of the eigenvector analysis. Fig. 13a shows 
the root locus out of ground effect. According to Ref. 1, at low speed the analysis 
yields two unstable modes which are the phugoid (Ph) and the lateral oscillation 
(L.O.), two aperiodic and stable modes such as the roll (R.S.) and pitch (P.S.) 
subsidence, and spiral (Sp) and heave (Hv) modes, that are both stable. Because 
of the vehicle axial symmetry about zb, the two modes Ph and L.O. degenerate 
in a single conical mode whose eigenvalues collapse as V — > 0. As the flight speed 
increases, the phugoid eigenvalue diminishes its imaginary part until to intersect 
the horizontal axis at about 8 m s -1 , where two unstable longitudinal aperiodic 
modes, i.e. A and B (see Fig. 13a), appear. As for L.O. and R.S., the corresponding 
eigenvalues are approximatively constant until to a velocity of about 20 m s _1 , 
whereas the pitch subsidence P.S. eigenvalue decreases, reaching its minimum at 
about 20 m s -1 . It is worth to remark that, in the eigenvectors calculated out of 
ground effect the Euler angles lag the other state variables and this is in agreement 
to the classical flight mechanics, where the Euler angles do not appear in the moment 
equations of motion. In ground effect the root locus presents some modifications 
which become more pronounced as h/R goes to zero. In particular, now, for h/R — > 
0, the eigenvectors associated to Sp, P.S. and R.S. show Euler angles that do not 
lag the other state variables. This is a peculiarity of all the vehicles that fly in 
proximity of the ground, which is caused by the fact that the R.H.S. of the moment 
motion equations depend on the Euler angles through the aerodynamic coefficients 27 
and thrust terms. For an assigned velocity, a significant reduction of the eigenvalues 
real part is observed for h/R — > 0. Because of the hull and rotors moments, under 
h/R < 2 Hv becomes an oscillating and stable mode whose eigenvalue imaginary 
part increases when h/R goes to zero, whereas the phugoid degenerates in a sort 
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of short period mode whose stability margin increases when the vehicle approaches 
to the ground. As for the lateral oscillation, it increases its eigenvalue imaginary 
part, and results to be stable for h/R < 2. According to the dynamics of the ground 
effect machines 28 ' 29 , the eigenvalues-eigenvectors analysis shows that at relatively 
low velocities, for h/R xs 0.5, all the modes result to be stable (See Fig. 13 d), due 
to the fact that the derivatives C m9 , Cm-po^Ci^ and C\ Tif) are negative. Therefore, 
the Euler angles play a role of paramount importance in the rotorctaft stability in 
ground effect, especially at low h/R. 

Conclusion 

This study analyzes the performance and stability of a shrouded-fan UAV by 
means of two models which determine force and moment developed by the vehicle 
in the different flight conditions. The models calculate force and moment developed 
by the hull and by the rotors and take into account the reciprocal interaction between 
fuselage and rotors and the ground effect. These features make the main advantage 
of the proposed method with respect to the other models known in the literature. 
The limitation of the method is that both the models are based on the potential 
flow theory where the flow structure about the airframe is in advance assigned and 
corresponds to the normal working state. The present study shows that the two 
models give results in good agreement with the various source of data existing in 
the literature. 

The trim analysis in and out of ground effect was accomplished and the corre- 
sponding results are discussed. Out of ground effect the same characteristics just 
known in the literature are recovered, whereas in ground effect sizable changing in 
the controls that are caused by the interaction between ground and rotorcraft are 
observed. The influence of the ground inclination with respect to the horizon is also 
studied which is responsible of seeable changing in the controls especially for what 



26 



concerns the lateral cyclic pitch. 

Finally the vehicle stability is analyzed through the eigenvalues-eigenvectors anal- 
ysis of the linearized motion equations. While the stability out of ground effect is 
in somewhat comparable with the data known in the literature, the hull and rotors 
force and moment, developed in ground effect, cause great modifications in the ve- 
hicle stability. According to the classical flight mechanics, the modes calculated out 
of ground effect exhibit Euler angles which lag the other state variables since such 
angles do not appear in the R.H.S. of the motion equations, whereas in ground effect, 
due to the incremental force and moment which depend on the Euler angles, the 
diverse eigenvectors show the Euler angles that do not lag the other state variables. 
This property, which yields stable modes in ground proximity, is in agreement with 
the behavior observed in the ground effect machines. 
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Table I. Characteristics of the UAV. 
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List of captions 

Fig. 1 Reference rotorcraft. 

Fig. 2 Rotor working regimes. 

Fig. 3 Vehicle aerodynamics in normal working state out of ground effect. 

Fig. 4 Vehicle aerodynamics in ground effect. 

Fig. 5 Aerodynamic coefficients out of ground effect. 

Fig. 6 Aerodynamic coefficients out of ground effect: 
Continuous lines for the present data. The symbols are from Ref. 23. 

Fig. 7 Aerodynamic coefficients in ground effect. 

Fig. 8 Thrust and torque coefficients in axial flight out of ground effect: 
dashed and continuous lines are, respectively, for free and shrouded rotors. 
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Fig. 9 Rotor in ground effect. 

Fig. 10 Influence of height, attitude and sink rate 
on the rotors characteristics. 

Fig. 11 Trim calculation in horizontal flight at different h/R. 

Fig. 12 Trim controls at h/R =1 with inclined ground wall. 

Fig. 13 Root Locus in horizontal flight at various h/R. 
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Figure 1: Reference rotorcraft. 
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Figure 2: Rotor working regimes. 
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Figure 3: Vehicle aerodynamics in normal working state out of ground effect. 
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Figure 4: Vehicle aerodynamics in ground effect. 
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Figure 5: Aerodynamic coefficients out of ground effect. 
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Figure 6: Aerodynamic coefficients out of ground effect: Continuous lines for the 
present data. The symbols are from Ref. 23. 
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Figure 7: Aerodynamic coefficients in ground effect. 
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Figure 8: Thrust and torque coefficients in axial flight out of ground effect: dashed 
and continuous lines are, respectively, for free and shrouded rotors. 




41 



Figure 9: Rotor in ground effect. 
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0.032 - Figure 10: Influence of height, attitude and sink rate on the rotors characteristics. 
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calculation in horizontal flight at different h/R. 
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Figure 12: Trim controls at h/R =1 with inclined ground wall. 




Figure 13: Root Locus in horizontal flight at various h/R. 
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